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Context. The search for B-mode polarization fluctuations in the Cosmic Microwave Background is one of the main challenges of 

modern cosmology. The expected level of the B-mode signal is very low and therefore requires the development of highly sensitive 

instruments with low systematic errors. An appealing possibility is bolometric interferometry. 

Aims. We compare in this article the sensitivity on the CMB angular power spectrum achieved with direct imaging, heterodyne and 

bolometric interferometry. 

Methods. Using a simple power spectrum estimator, we calculate its variance leading to the counterpart for bolometric interferometry 

of the well known Knox formula for direct imaging. 

Results. We find that bolometric interferometry is less sensitive than direct imaging. However, as expected, it is finally more sensitive 

than heterodyne interferometry due to the low noise of the bolometers. It therefore appears as an alternative to direct imagers with 

different and possibly lower systematic errors, mainly due to the absence of an optical setup in front of the horns. 
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Introduction 

The detection of primordial gravity waves through B-mode po- 
larization anisotropies in the Cosmic Microwave Backgroud is 
one of the most exciting challenges of modern cosmology. It 
could provide direct information on the energy scale of infla- 
tion, possibly associated with GUT (Liddle and Lyt h, 200TJ| . It 
would also allow one to investigate the standard cosmological 
model in detail through consistency tests involving the spectral 
indices of scalar and tensor perturbations and their amplitude ra- 
tio (|LangIoIs720M|. 

Despite the weakness of the expected signal, many teams 
have decided to join the quest for the B-modes and to construct 
dedicated instruments that must combine exquisite sensitivity 
and precise control of systematic effects. Most of the projects 
proposed up to now use direct imagers, a concept that has proven 
to be very sensitive. However they might be affected by signif- 
icant systematic effects such as ground-pickup and beam differ- 
ences that would less affect an interferometer having no optics 
before the entry horns. Thus, we investigate the possibility of de- 
veloping a high sensitivity interferometer dedicated to B-mode 
searches. A bolometric interferometer would combine the high 
sensitivity of bolometers with the clean optics of an interferom- 
eter and could complement the ongoing imaging projects. 

In this article we investigate the sensitivity achieved by such 
an instrument and compare it with direct imagers and heterodyne 
interferometers. In section [T] we review useful quantities regard- 
ing interferometry in general and the reconstruction of visibili- 
ties with a bolometric interferometer. We define in section [2] a 
simple power spectrum estimator under the assumption that E 
and B visibilities can be extracted from Q and U Stokes param- 
eter visibilities (this is in itself an important issue and is beyond 



the scope of this article). From this estimator we derive a formula 
for the C( uncertainty with a bolometric interferometer equiva- 
lent to the famous one derived in (Knox, 1997) for imagers. We 
compare the sensitivities achieved with imagers and bolometric 
and heterodyne interferometers in section [3] We discuss our re- 
sults in section|4] 



1. Definitions and assumptions 

1.1. Useful interferometric quantities 

Here, we discuss basic definitions regarding interferometry in 
general and bolometric interferometry in particular. An interfer- 
ometer observes the visibilities of an incoming radiation field 
S (n) that are defined as: 



v«(») = I 



S (n)A(n) exp(2/7r u ■ n)dn 



(1) 



where u is the baseline defined as the vector separation between 
the horns D in units of the electromagnetic wavelength of the ra- 
diation u = D/A. The way visibilities are actually reconstructed 
in a bolometric interferometer is described in section 11.21 The 
beam of the input horns A(n) is normalized to one at maximum. 
In the flat-sky approximation, one can write the visibility as a 
convolution of the Fourier modes of the incoming radiation with 
the Fourier transform of the input beam: 



v s (w) = I 



§ (u)A(u - v)dv = S («) <8 A(«) 



(2) 
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where the " denotes the Fourier transform. Note that outside the 
flat-sky approximation, a similar formula is to be expected, al- 
though more complicated (B unn and W hite, 2007). This convo- 
lution expression states that the sky cut in real space performed 
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by the input beam is expressed as a convolution in Fourier space. 
The correspondance with multipoles is 6 — 2im. Each horn cov- 
ers a solid angle Q defined as: 



Q 



/ 



A(n)dn. 



(3) 



If we approximate the beam by a Gaussian defined by its RMS 
cr, its FWHM is cr x 2 V21n2 a* 2.35cr and in that case the solid 
angle subtended by the horn is Q. - Ino 2 which can be related 
to the fraction of the sky observed with the horn f^y - Q./4-n. 
We assume that the horns are placed on a square grid so that 
the baselines are also located on a square grid. The minimum 
spacing between two horns is obtained by packing at a distance 
equal to their diameter. As they are diffraction limited, their sec- 
tion, solid angle and wavelength are such tha{H5'Q =* A 2 so that 
their distance is D/, = 2/1/ y/nCl. The spacing between visibilities 
in Fourier space is therefore: 



1 



(4) 



scheme is such that equivalent baseline^jare summed coherently 
- they correspond to the same phase difference. In this case, the 
noise covariance matrix on the complex reconstructed visibilities 
is diagonal and has the form: 



Nu 



St) 



4 NET 2 Q 2 JV/, 1 



N, 



a£j(0' 



(8) 



The first factor 2 comes from the fact that when measuring polar- 
ized visibilities, one cannot have access to both Q and U Stokes 
parameters at the same time (Charlass ier et al., 2 008). The sec- 
ond is due to the fact that here we are dealing with the covariance 
matrix of the complex visibilities instead of their real or imag- 
inary parts as in (Charlassier et al., 2008). Nf, is the number of 
entry horns, N, is the number of time samples, N eq (i) is the num- 
ber of equivalent baselines corresponding to baseline / and NET 
is the noise equivalent temperature expressed in fiK/ y/Tlz. An 
extensive analytical and Monte-Carlo based study of the recon- 
struction of the visibilities in bolometric interferometry can be 
found in (Charlassie r et al., 2008| l. 



In our Gaussian approximation, the Fourier transform of the pri- 
mary beam is: 



A(u) = Qexp(-7tQH 2 ) 



(5) 



so that in terms of a Gaussian in Fourier space, the resolution is: 



1 



V2~^Q 



(6) 



The baselines are separated by u m ; n and each point in baseline 
space has a resolution <x u , therefore: 



2V2^3. 



0"u 



(7) 



We can conclude that the different baselines are almost inde- 
pendent. In reality, the distribution of the electric field vanishes 
outside the horn aperture so that the primary beam is not an ex- 
act Gaussian and the beam in Fourier space is also truncated at a 
radius of u m i n (Whi te et al., 19 99). This does not change the fact 
that different baselines are almost independent. 

1.2. Visibility reconstruction in bolometric interferometry 

In a heterodyne interferometer, the visibilities are directly ob- 
tained using a correlator that gives as an output the correlated 
signal E\E* coming from two antennae in a coherent way. 
Bolometers are incoherent detectors that measure the time av- 
eraged incoming power. Bolometric interferometers are there- 
fore additive interferometers where the visibilities are obtained 
from squaring the sum of the signal coming from two horns : 
P = (|£, + E 2 \ 2 ) = (|£!| 2 ) + (|£ 2 | 2 ) + 2(E l E*). When a large 
number of horns are used, the signal detected by the bolometers 
is a linear combination of all available visibilities. The use of 
phase shifters on each of the input channels allows one to recon- 
struct the complex visibilities of the four Stokes parameters V/, 
Vg, Vj/ and Vy We have shown (Charlass ier et al., 2008| l that 
this reconstruction is done optimally when the phase-shifting 



1 Actually, for a Gaussian beam, the illumination function on the en- 
try of the horn is not flat (but is Gaussian) and the size of the horn has 
to be larger that what is quoted here by a factor of about 2. 



2. Power spectrum estimator 

2.1. E and B fields from the Stokes parameters 

In the flat-sky approximation, the E and B polarization fields are 
related to the Stokes parameters by a simple rotation of angle 0, 
the angle between u and the u x axis (Zaldarriaga, 2001). In terms 
of visibilities, this can be written: 



Vq(m) = J cos 2<p E(v) - sin 2cf> B(v) 
\/ u (u) = J sin 20 E(v) + cos 20 B(y) 



A(u - v)dv 
A(u - v)dv 



(9) 



For the simplicity, we do not discuss the E/B separation here and 
assume that one can obtain a set of pure E and B visibilities from 
the Stokes parameter visibilities defined aJM 



V £ (m) = J E(v)A(u - v)dv 
y B (u) = f B(v)A(u - v)dv 



(10) 



and the covariance matrix of the B visibilities contains the BB 
angular power spectrum: 



(Vb(m)V*(h')) = f(B(v)B*(v'))A(u-v)A*(u' -v')dvdv(ll) 

= fc BB (v)A(u-v)A*(u'-v)dv (12) 

= 6(u-u')x fcf(v) |A(H-v)| 2 dv (13) 

the last equation comes from the fact that the different baselines 
we measure are independent from the beam point of view. 

In the presence of noise and assuming the power to be flat 
enough to be taken out of the integral (recall that I - 27ru): 

(v b (m)V*(h')) = S(u - u) x Cf B f |i(v)| 2 dv +N(u, u). (14) 

=n/2 

2 Equivalent baselines are sometimes called redundant baselines in 
the literature. They correspond to different pairs of horns separated by 
the same vector and therefore corresponding to the same point in (u, v) 
space. 

3 Obtaining such pure E and B modes is a complex issue in itself and 
deserves a full study. 
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As said before, the noise covariance matrix of the Stokes param- 
eter visibilities is diagonal and we assume that it is still the case 
for that of Vg and Vg. Labeling u and u' with indices i and j and 
Vb(h) and Vb(h') by V, and Vj, one gets: 



(v,v;) 



Q. 

Cc-Sij + Nijdij. 



(15) 



2.2. A simple pseudo-power spectrum estimator and its 
variance 

From the above equation, it is obvious that the simplest unbiased 
estimator of the power spectrum is: 



C,= 



2 1 



J] (V,Vf-M,) 



(16) 



where N±(€) is the number of different baselines corresponding 
to multipole I. The variance of this estimator is (Cc is the true 
power spectrum): 



Var(C f ) = (C 2 ( ) - C] 



W*(€) 



^(V,V* - Nid 



-C 



N: 



k?( C 



' + l Ni - 



(17) 
(18) 

(19) 



where we used Wick's theorem to calculate the fourth order mo- 
ments and the fact that each of the N+(€) different baselines 
contributing to Q is measured independantly so that their vari- 
ances add linearly. If one makes the additional assumption that 
all of these different baselines have the same noise variance 
Nij - <Ty 6ij, the error on the power spectrum reads: 



AC, 



1 



'N+(fi 



Cc 



2<t 2 \ 



(20) 



which is the equivalent for interferometry of the well known 
imaging-oriented Knox formula (Kn ox, 1997] >. The expression 
for N±{() is the number of different modes one can have access 
to at a given I, We assume that we are considering a bin in vis- 
ibility space Au = M/ln centered at u; the number of modes is 
the ratio between the available surface of the bin miAu (we only 
consider the top part of the Fourier plane as the modes in the 
bottom part are the same) to the effective surface of the beam in 
Fourier space 27rcr^ : 



LTicri 



(21) 



we therefore find the same formula as for an imager, except for 
the noise part of course: 



AC? l = 



' 2iU y M 



C ( + 



2 < 



(22) 



where the noise on the visibilities is taken 
from dCharlassier et al, 2008| l. 



_BI 




(23) 



3. Comparison with an imager and a heterodyne 
interferometer 

3.1. Analytical formulae 

The above expression is the same for both heterodyne interfer- 
ometry and direct imaging, only the expression of <x v changes. 
For heterodyne interferometry, if the noise equivalent tempera- 
ture of one of the two input channels of the correlator is NEThi, 
the noise on the reconstructed Stokes parameter visibility cal- 
culated with N, time samples and averaged over N eq equivalent 
baselines is given by: 



2 V2 NEThi O 



(24) 



The first factor 2 comes from the multiplication of the two 
sine waves, a factor V2 from the fact that two correlators 
are involved when calculating Stokes parameters visibilities, 
another factor V2 appears because we are talking about the 
noise on the complex visibility instead of its real or imaginary 
part. Finally a factor 1 / V2 is regained because two sets of in- 
dependent measurements of the Stokes parameters visibilities 
can be simultaneously obtained if one forms all the possible 
complex correlations. The expression we find is in agreement 
with ( jHobson a nd Magueijo, 1996} and dWhite et al., 1999| >. In 
the direct imaging case, the error on the power spectrum is taken 
from (Kn ox, 1997| l and adapted to partial sky polarized measure- 
ments: 



AC™ = 



'(2{+l)f sky M 



C f + 



4NETj 2 m Q^ 

N h B 2 c N, 



(25) 



where Bg - exp(-£ of /2) is the imager's beam transfer func- 
tion. In the imaging case Q is of course defined as the solid an- 
gle covered on the sky Q = 4-nf^y (the fact that the integral of 
the primary beam is the total solid angle covered on the sky is 
specific to interferometry). Note that the factor of 4 is obtained 
by a factor of 2 on the polarized NET for polarization sensitive 
bolometers and another factor of 2 due, as before, to the fact that 
Q and U cannot be obtained at the same time. 

As the sample variance term is exactly the same whatever 
technique is used (as expected), we are only interested in com- 
paring the noise terms. We assume in the following that we are 
comparing three instruments observing the same fraction of the 
sky / s ky from the ground for the same duration: 

- a direct imager with N/, horns, an angular resolution given 
by cr beam and a NET lm ^150 /jK/ VHz, as stated for Clover 
in dNorfh et al, 20071 ) at 97GHz in Chile. 

- a heterodyne interferometer with a primary beam cov- 
ering / s ]<y, using a square array of N/, input channels 

ea ch with NET m ^ 250 //K/ Vllz, as stated for QUIET 
in dSamtleben et al, 2007) 1 at 90 GHz in Chile. 

- a bolometric interferometer with the same characteristics as 
the heterodyne one but with a NET identical to that of a bolo- 
metric imager NET BI ^150 /iK/ VHz. 

The choice of 90GHz is motivated by the fact that a packed array 
of 20 degrees FWHM primary horns interferometer operating at 
these frequencies would cover the multipole range relevant for 
primordial B-mode signals (25 < I < 200). It is also at these 
frequencies that coherent and bolometric detectors can operate 
simultaneously. We will use the direct imager as a reference and 
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30 



and actual val- 



Figure 1. Approximate formula (red line) from Eq. 

ues (black dots) for the number of equivalent baselines as a function of 

I = 2im for 256 entry horns on a square compact grid. 



calculate the ratio of the direct imager's noise error to that of the 

interferometers. This ratio therefore should be less than one if 

the direct imager is more sensitive from the strict noise point of 

view. 

For the bolometric interferometer, one gets: 



AC 



i m 



ACf 



1 (N eq Y 1 / NET lm 
2 X UJ X ^ X lNET, 



2 \N h 



1 
3 



BI 



(26) 
(27) 



as the NET are the same for the bolometers used for imaging or 

for bolometric interferometry. 

For heterodyne interferometry, one gets: 



AC ' m 



AC^ 



HI 



1 



N, 



eq 



N, 



1 /NET 



X B? X \NET 



i in 



HI 



(28) 



We can see that the ratio N eq /Nh is always less than one, which 
gives a clear advantage to direct imaging from the strict point 
of view of the noise. This ratio appears squared in the ratio of 
imaging to bolometric interferferometry and without power in 
the ratio of imaging to heterodyne interferometry but in the latter 
case, the NET ratio is also less than one, penalising heterodyne 
interferometry. 

3.2. Approximate expressions and simulations 
The number of equivalent baselines for a square horn array is: 

u. 



N eq = 



N h - 



N h - 



(29) 



If one averages over directions in the baseline plane at a given 
\u\, a good approximation of N eq as a function of I is given by 
(see Fig. [Til: 



Aeq ^ x V2 

N h ~ 2VAa 



where 



2VAa 



(30) 



One finally finds that a good approximation of the sensitivity 
ratio is: 



ACj m 



ACf 



BI 



- 1- 



V2 



2^ 



1 

2~Bj 



(31) 



and: 

AC 1 /" 



ACf 



* 1- 



V2 

2VA^ 



n 



'«/ 



1 / NET lm 
X 4B? X lNET HI 



(32) 



These approximate formulae have been compared with actual 
calculations of the number of equivalent baselines for square 
arrays. We have chosen 256 horns for the comparison and we 
compare bolometric and heterodyne interferometers with im- 
agers having a low angular resolution of one degree, BICEP- 
like ( Yoon et al, 2006) and a high one of 10 arcminutes, Clover- 
like (North et al, 2008). The results are shown in Fig. [2] We have 
chosen to only consider the multipole region between and 200 
as for higher multipoles, interferometers are less sensitive due 
the loss of coherence between largely separated horns. Note that 
the effect of coherence loss for the long baselines and the band- 
width smearing have not been taken into account here and might 
have a significant effect. 

4. Discussion 

The sensitivities of the three different techniques only differ in 
the way the instrument filters the multipoles observed in the sky. 
An imager is affected by its resolution on the sky while an inter- 
ferometer is affected by the ratio between the number of equiv- 
alent baselines and the number of horns as a function of multi- 
poles. All of these filtering factors are less than one. However, 
imagers are usually operated in such a way that they are not lim- 
ited by their angular resolution in the multipole region of inter- 
est, in that case Be =* 1, and the imager is always more sensi- 
tive than an interferometer (bolometric or heterodyne). From the 
strict point of view of sensitivity, interferometers can therefore 
only compete with low angular resolution imagers. 

There is a large difference in sensitivity between bolomet- 
ric and heterodyne interferometers compared to an imager: the 
ratio Neq/N/, acts quadratically on the variance for a bolomet- 
ric interferometer while it acts linearly for a heterodyne instru- 
ment. This is due to the fact that with a heterodyne interferom- 
eter, equivalent baselines are averaged after their measurement, 
resulting in a 1/A^q factor on the variances. In a bolometric in- 
terferometer, the signals from all Nh horns are added together 
multiplying the noise variance by Nh while the coherent sum- 
mation of equivalent baselines performs an efficient l/N'L re- 
duction of the noise. This finally results in a factor Nh/2N eq for 
the variance of a bolometric interferometer relative to a hetero- 
dyne one. This is largely compensated by the difference in NET 
between bolometric instruments and coherent ones. When com- 
paring them, the ratio of their NET also appears quadratically 
and favours bolometric instruments that are dominated by the 
photon noise rather than by that of the amplifiers. This situation 
may change in the future with the improvements of the HEMT 
technologies but at frequencies around and above 100 GHz we 
are unlikely to face photon noise limited HEMTs in the near fu- 
ture. The difference between the NET would be even greater in 
space where the bolometers NET would drop as the background 
temperature while that of the coherent instruments would remain 
roughly constant. 

With the present technologies of bolometers and coherent 
amplifiers, the hierarchy in terms of sensitivity between the three 
techniques (and layout) studied here is very clear for the multi- 
pole range 25 < ( < 200 where the primordial B-mode signal 
is expected to be maximal. Imagers are the most sensitive, bolo- 
metric interferometers have a lower sensitivity, the ratio drop- 
ping quadratically with the multipole considered. Heterodyne 
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N h =256 - FWHM lm = 0.17 - FWHM hF 20.0 - FWHM bF 20.0 



N h =256 ■ FWHM lm = 1.00 ■ FWHM hF 20.0 ■ FWHM BI = 20.0 



1.2 



Direct Imager I Heterodyne Interferometer 

Direct Imager I Bolometric Interferometer 

Noise only 

Sample variance included 




100 
I multipole 



1.2 



1.0 



a 0.8 



1 



0.6 



0.4 



0.2 



0.0 



Direct Imager I Heterodyne Interferometer 

Direct Imager I Bolometric Interferometer 

Noise only 

Sample variance included 



NET lm = 150nK/Hz" : 
NET m =150\iK/Hz" ! 
NET m = 250\iK/Hz K 
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Figure 2. Power spectrum error bar ratio between an imager and a bolometric interferometer (red solid) and between an imager and a heterodyne 
interferometer (blue dashed). The thick lines only include the noise contribution while the thin lines also include sample variance (with a tensor to 
scalar ratio of 0. 1 ). In each case we have assumed 256 horns, a sky fraction defined by the 20 degrees FWHM primary beam of the interferometers. 
The left panel corresponds to a 10 arcminute resolution imager and the right one to a one degree imager. The points are obtained with a precise 
measurement of the number of equivalent baselines and the lines are from the approximate expression given above. The green dotted line shows the 
region where the relative sample variance AC?/C{ = J2/(2tf !iky A{) is greater than one. The ratio of the heterodyne interferometer to the imager 
depends on the square of the NET ratio; we have assumed 250 and 150 mK/ VHz taken from ( Samtleben et al, 2007 ) and Clover ( North et al, 2007 ) 
at ~90 GHz for Chile. It is straightforward to scale the curves for different values of the NET. 



interferometers have an even lower sensitivity but the ratio with 
an imager drops less rapidly. They remain however less sensi- 
tive than bolometric interferometers in the range of multipoles 
considered here, where the largest primordial B-mode signal is 
expected and where the lensing of the E-modes into B-modes 
is still low. At higher multipoles however, the bandwidth smear- 
ing effect and loss of coherence would be a real issue for bolo- 
metric interferometers while in a heterodyne interferometer, the 
separation into small bands would prevent the sensitivity from 
dropping. 

The main remaining question is whether the gain in terms 
of systematic effects is worth the price of this sensitivity re- 
duction if one builds an interferometer instead of an imager. In 
terms of optics for instance, an interferometer directly observes 
the sky. The primary beam is therefore only set by that of the 
horns, while in an imager, the telescope (mirror or lenses) pro- 
duces sidelobes inducing poorly predictible ground pickup that 
often prevent one from reaching the nominal sensitivity. An in- 
terferometer is also completely insensitive to spatially uniform 
polarized signals that vary with time such as polarized atmo- 
spheric contamination. These could also prevent an imager from 
reaching its nominal sensitivity by adding some spread in the 
noise. These examples mitigate the statistical sensitivity loss 
of an interferometer with respect to an imager. The differences 
in terms of systematic effects between imagers and bolomet- 
ric and heterodyne interferometers are not obvious and deserve 
a detailed quantitative study in continuation of the work done 
by dBunn, 2007] l. 
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